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We discuss a new class of solutions to the Einstein equations which describe a primordial black 
hole (PBH) in a flat Friedmann background. Such solutions arise if a Schwarzschild black hole is 
patched onto a Friedmann background via a transition region. They are possible providing the black 
hole event horizon is larger than the cosmological apparent horizon. Such solutions have a number of 
strange features. In particular, one has to define the black hole and cosmological horizons carefully 
and one then finds that the mass contained within the black hole event horizon decreases when 
it is larger than the Friedmann cosmological apparent horizon, although its area always increases. 
These solutions involve two distinct future null infinities and are interpreted as the conversion of 
a white hole into a black hole. Although such solutions may not form from gravitational collapse 
in the same way as standard PBHs, there is nothing unphysical about them, since all energy and 
causality conditions are satisfied. Their conformal diagram is a natural amalgamation of the Kruskal 
diagram for the extended Schwarzschild solution and the conformal diagram for a black hole in a 
flat Friedmann background. In this paper, such solutions are obtained numerically for a spherically 
symmetric universe containing a massless scalar field, but it is likely that they exist for more general 
matter fields and less symmetric systems. 

PACS numbers: 04.70.Bw, 97.60.Lf, 04.25.Dm, 95.35.+d 

I. INTRODUCTION 

In a recent paper (henceforth Paper I) we studied numerically the growth of primordial black holes (PBHs) in 
a universe containing a massless scalar field but no other matter. Following Hamade and Stewart 0], the double- 
null formulation of the Einstein equations was used. This is a powerful tool for investigating regions outside the 
cosmological horizon and inside the black hole horizon simultaneously. On the assumption that the PBH is formed 
from a local initial density perturbation which propagates causally, the black hole was modelled by matching a 
Schwarzschild solution to an exact flat Friedmann solution across a null surface. Initial data were specified on an 
outgoing and ingoing null surface. They were assumed to be exactly Friedmann on the outgoing surface and outside 
the matching boundary on the ingoing surface but some perturbation of Friedmann inside the matching boundary. 

In all the solutions considered in Paper I, the black hole event horizon (BHEH) was assumed to be smaller than the 
cosmological apparent horizon. However, in some circumstances (e.g. in the inflationary scenario), the perturbation 
- and perhaps even the black hole itself - may extend well beyond it. In this case, the only upper limit on its size 
comes from requiring that the perturbed region be part of our universe rather than a separate closed universe. The 
condition for this has been derived precisely for the situation in which the collapsing region is homogeneous and the 
equation of state is p = fcp 0. In particular, this includes the massless scalar field case, since this is equivalent to a 
stiff fiuid with k = 1 providing the gradient of the field is everywhere timelike and vorticity free (as expected). The 
peculiarity of the formation of PBHs in a universe with a stiff fiuid was discussed in Ref. [J , and this suggests that 
the gradient of the scalar field could have an important hydrodynamical effect. 

We study such "super-horizon" solutions in this paper, although our numerical results only cover the scalar field 
case. It turns out that these scalar field solutions have some rather strange properties. For example, the mass 
contained within the BHEH decreases when the PBH is larger than the cosmological apparent horizon, although 
the area always increases. Since our system satisfies all energy conditions, the mass decrease of super-horizon PBHs 
we show here is completely different from the black hole mass decrease due to phantom energy accretion Also 
the conformal diagram for these solutions is interesting, being a natural extension of the Kruskal diagram to the 
cosmological context. The scalar field usually has a timelike gradient in the numerical simulations shown here and 
this means that it is equivalent to a stiff fliud. It is likely that there are analogues of these solutions for more general 
fluids with p — kp providing k > 1/3. However, we do not study these solutions here. In all these cases, we need to 
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define cosmological and black hole horizons very carefully when their sizes are all of the same order. 



II. FORMULATION 



A. Double-Null Formulation of Einstein equations 



As in paper I, we consider a massless scalar field in general relativity, for which the stress-energy tensor is 

Tab - - isab*"*^. (2.1) 

The Einstein equations are 

Rab - \9abR = &7^Tab, (2.2) 

and the equation of motion for the scalar field is 

□ * = = 0. (2.3) 

We focus on a spherically symmetric system, for which the line element can be written in the form 

ds^ = -a^{u, v)dudv + r^{u, v){de^ + sin^ ed<f)^), (2.4) 

where u and v are advanced and retarded time coordinates, respectively, a is the metric function and r is the "area 
radius" (the proper area of the sphere of constant r being Airr'^). Equations H2.2|l and H2.3|l then imply that we have 14 
first-order partial differential equations and two auxiliary equations. These equations are given explicitly in Section 
2 of 0- We adopt units in which G = c = 1. 



B. Misner-Sharp mass and trapping horizons 

The existence and position of marginal surfaces can be inferred from the form of the Misner-Sharp mass This 
is a well-behaved quasi- local mass in spherically symmetric spacetimes |2| , which can be written as 

,„.L(i + '-^)^ (2^5) 
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Combining this with the equations in Section 2 of 0], we can derive the following useful relations: 

8^r2r„(*,„)2 
= 5 , 

8^r2r„(4',,)2 



(2.6) 

(2.7) 
«- 

In this paper, we adopt the "trapping horizon" framework introduced by Hayward 0, |^ because this provides a 
systematic and mathematically transparent view of the sort of cosmological black holes treated here. In the context 
of this paper, trapping horizons and conventional apparent horizons are almost equivalent but they need not be in 
more general situations. 

Although the values for r^y and are not geometrical invariants, their signs are, which leads to the following 
definitions, (i) A metric sphere is said to be trapped if r^^r^u > 0. A sphere with r^u < and r^^ < is future trapped, 
while one with r^u > and r^y > is past trapped, (ii) A metric sphere is said to be untrapped if r^yV^^ < 0. On an 
untrapped surface, dy is outgoing if r^y > and ingoing if r^y < 0. More generally, a spacelike or null normal vector 

is outgoing if z'^r^a > and ingoing if z°T^a < 0. (iii) A metric sphere is said to be marginal if r^yV^u = 0. A sphere 
with r^y = is future marginal if r,„ < 0, past marginal if > and bifurcating marginal if r,„ — 0. It is described 
as outer marginal if r^uv < 0, inner marginal if r^^v > and degenerate marginal if r^uv = 0. It is easily seen that a 
sphere is marginal if and only if r = 2m, trapped if and only if r < 2m and untrapped if and only if r > 2m. 

The closure of a hypersurface foliated by a future or past, outer or inner marginal sphere is called a (nondegenerate) 
trapping horizon. In Hayward's approach, the black hole apparent horizon is replaced with a "future outer trapping 
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horizon" (FOTH), while the cosmological (white hole) apparent horizon is replaced with a "past outer trapping 
horizon" (POTH). There is a critical difference between an apparent horizon and a trapping horizon. An apparent 
horizon is defined only on a prescribed spacelike hypersurface, which is usually required to be a Cauchy surface, so it 
exists only if the spacetime is strongly asymptotically predictable On the other hand, a trapping horizon is the 
trajectory of a marginal surface in the whole spacetime, and it is defined locally whenever a null foliation is possible. 
The spacetime does not need to be strongly asymptotically predictable. It should be stressed that Hawking's "area 
theorem" 1^ only refers to the BHEH. On the other hand, Hayward has shown that it also applies for a FOTH 
providing the null energy condition holds. 

Providing the black hole horizon is within the cosmological horizon and u and v are the standard double-null 
coordinates in the asymptotic Friedmann region, the FOTH and the POTH correspond to the conditions = and 
r K = 0, respectively. However, the situation is more complicated if the black hole horizon is outside the cosmological 
horizon. In this case, we can still define trapping horizons by the conditions ~ and = 0, but these are 
no longer everywhere identified with a FOTH and a POTH. Using the above equations, we can show that along a 
trapping horizon, on which r = 2m, we have 

[aV,„ + 167rrV,„(*^„)2] du + [aV,^, + 167rrV,„(«',„)^] dv = 0. (2.8) 

On trapping horizons, which have r ^, = and r.„ = 0, we therefore have 

a^du + 16TTr^{^^y)^dv = 0, (2.9) 
16TTr'^{-i'^ufdu + a'^dv = 0, (2.10) 

respectively. We conclude that trapping horizons are non-timelike in this system. More precisely, a trapping horizon 
with r^y = has u = const if and only if y = 0, while a trapping horizon with r^y = has v = const if and only if 
^^11 = 0. Except for these special cases, trapping horizons are spacelike. 

The form of a{u, v) and r(u, v) in the exact flat Friedmann model is given in Appendix B of Paper I. This implies 
that the cosmological particle horizon has u — 0, while the cosmological apparent horizon, which is a POTH, has 
3u + V — 0. This shows that the cosmological apparent horizon is spacelike and outside the particle horizon. The 
conformal diagram of the spacetime is indicated in Fig. ^ which shows the initial (big bang) spacelike singularity. The 
cosmological apparent horizon also coincides with the Hubble horizon in this case. If one has a black hole embedded 
in an exact or asymptotically flat Friedmann model and smaller than the cosmological particle horizon, the conformal 
diagram will change to the form indicated in Fig. 13 This now contains a BHEH, a FOTH and a final (black hole) 
spacelike singularity. 

III. INITIAL DATA FOR PBHS 

A. Structure of initial data 

The initial data are prescribed on the two null surfaces u = uq and v = vq, with the region of calculation being the 
diamond [uo,ui] x [uo,fi]. We have three independent functions on the two null surfaces: a, 5* and Two of them 
can be chosen freely and the third is determined by the constraint equations on the null surfaces. It is convenient to 
choose 

a(uo, w), a(u, vq), *(wo, «), '^{u, vq) (3.1) 

as the free initial data and to regard 

r{uQ,v),r{u,vo) (3.2) 

as being determined by the initial value equations. We can regard ^(uo,u) and \l/(it, wq) as the physical degrees of 
freedom in the initial data, while the choice for a{uQ, v) and a{u, vq) fixes the gauge. 

In the flat Friedmann region, we adopt the coordinate system given in Appendix Bl of Paper I and impose flat 
Friedmann initial data for a and 5* on m = wq for vq < v < vi. On ?; = uq, we also choose flat Friedmann data for a 
for Mo < u < ui. For we use the same data on the initial null surface v = vq for uq < u < but ^' — const for 
Mm -l- Au < u < ui. This is equivalent to Schwarzschild data in coordinates penetrating the black hole. The sudden 
transition from flat Friedmann data to Schwarzschild data results in a discontinuity at u = which reduces the 
numerical accuracy. Hence we smooth the transition with some smoothing length Au; we use a quadratic function 
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between Um and Mm + Au, so that ^' and are continuous. More precisely, we impose the following initial data for 
a and s = V^^: 

(3.3) 

[u < Uni) 

{Uni < U < Uni + Au) , (3.4) 
{u > Mm + Am) 

on the initial null surface v = vq, and 

aHuo,v) = C'(^^y (3.5) 

s{uo,v) = ^In (^ ""2*"^ ^ +^^0, (3.6) 

on the initial null surface u = uq. Here C and sq are constants and, without loss of generality, we can choose C = 1 
and So = 0. Note that although one has a Schwarzschild vacuum for u > + Au, this situation only applies 
instantaneously at v = vq, since the inflowing matter will fill this up immediately. 

Figure |31 depicts the initial data for our numerical simulations, these being completely determined by the two 
parameters Mm and Am. When the matching region is very narrow, the BHEH is smaller than the Friedmann 
cosmological apparent horizon if Mm > mcah(i^o)- This corresponds to Fig. |2{a), where u — mcah(w) denotes the 
trajectory of the cosmological apparent horizon. However, if Mm < itCAH('yo)i the BHEH is larger than the Friedmann 
cosmological horizon, which corresponds to Fig. |3{b). Note that the cosmological apparent horizon would always go 
within the matching radius if one took vq small enough, so when the black hole forms is crucial. Note also that when 
we refer to a POTH in the background Friedmann model, we refer to it as the (Friedmann) cosmological apparent 
horizon. As discussed later, this relates to the distinction between a black hole that forms from collapse and an eternal 
black hole that exists "ab initio" . 



IV. RESULTS 

A. Black hole event horizon 

Our numerical code is based on that of Hamade and Stewart ^ with a modification, described in Appendix A 
of Paper I, to ensure greater accuracy. The initial data are prescribed on the two null surfaces v = vq = I and 
M = Mo = —2/3, this also fixing the units. The calculated region is the diamond contained by m = mo, m = t;o, m = mi 
and V = vi. On the initial null surface v = vq, we make the matching at m = Mm and use the smoothing length Am. 
As time proceeds, the m = const null rays will become ever more sensitive to r near the BHEH, so the calculation 
is stopped at f = Ml, when the null rays become too coarse to be resolved. The parameters used are summarised in 
Table ^ for six models. The BHEH is identified as the critical null ray m = mbheh, the null rays with m < mbheh 
going to r = oo and those with u > mbheh returning to r = 0. Therefore, the BHEH is only identified at the end 
of the calculation. Although the identification of the BHEH is rather imprecise, because of numerical errors and the 
finiteness of the calculated regions, it suffices for the physical interpretation of the results given below. It is found 
that all six models have a BHEH and Table indicates the initial ratio of the sizes of the BHEH and the Friedmann 
cosmological apparent horizon at m = mo- This complements Table I of Paper I, which only shows cases (A to D) in 
which the ratio is less than 1. All models except E have an initial BHEH larger than the Friedmann cosmological 
apparent horizon. In Model E the initial BHEH is slightly smaller than the Friedmann cosmological apparent horizon. 

B. Location of horizons 

Figure 21 shows the locations of the BHEH and trapping horizons in the {u,v) plane. It also indicates the signs of 
r,v and r,„. The region is future trapped if the signs are (— , — ), past trapped if they are (+, +) and untrapped if they 
are (+, — ) or (—,+). Trapping horizons may have either r.„ = or r = 0. 



a^(M, mq) = 



2 ( U + Vo 
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s(m. Mo) = < — 



V3 
2 

2 



(Am) — (Mm + Am — m) f + vq 



2AM(Mm + Am) 



In 



So 



Am 

2 (Mm + Mo) 



In 



So 



5 



TABLE I: Model parameters and the initial mass ratios of BHEH to the Friedmann cosmological apparent horizon. 



Models 


Ui-n 


Au 


Uo 


Ul 


vo 


Vi 


m-BHEn/mcAH 


E 


-1/2 


0.4 


-2/3 


2 


1 


4.5 


0.969 


F 


-1/2 


0.02 


-2/3 


2 


1 


4.5 


2.12 


G 


-1/2 


0.05 


-2/3 


2 


1 


4.5 


2.01 


H 


-1/2 


0.2 


-2/3 


2 


1 


4.5 


1.52 


I 


-0.6 


0.02 


-2/3 


2 


1 


4.5 


3.58 


J 


-0.4 


0.02 


-2/3 


2 


1 


4.5 


1.31 



There are two qualitatively different classes of models. In the first class, which includes all models except E, the 
initial BHEH is larger than the Friedmann cosmological apparent horizon. As seen from the figure, as v increases 
from vq = 1, two trapping horizons, one with r „ = and the other with r „ = 0, appear and cross each other in the 
future of the BHEH. In terms of the evolution with respect to v, before the crossing the trapping horizons with r = 
and r.„ = correspond to a POTH and a FOTH, respectively. After the crossing, the situation is reversed, so the 
trapping horizons with r^y = and r.„ = correspond to a FOTH and a POTH, respectively. After further evolution, 
the POTH crosses the BHEH. Thereafter it will coincide with the Friedmann cosmological apparent horizon, while 
the FOTH corresponds to the black hole apparent horizon. 

In the second class of models, which includes E, the initial BHEH is smaller than the Friedmann cosmological 
apparent horizon. Trapping horizons with r^y — and — appear and remain in the future and past of the BHEH, 
respectively. These two trapping horizons do not cross each other and always correspond to the FOTH and POTH, 
respectively. 



C. PBH mass change 

The area radii of the BHEH and trapping horizons for these models are shown in Fig. |S| It is seen that the area 
of the BHEH always increases for these models, which is consistent with the black hole area theorem If 
initial BHEH is larger than the Friedmann cosmological apparent horizon, i.e., for all models other than E, the areas 
of the trapping horizons with both = and r_„ = first decrease as v increases and then increase after crossing 
each other. It is interesting that, in terms of v, the two trapping horizons cross at a radius slightly larger than the 
BHEH at the crossing time. After the BHEH enters the future of the POTH, it soon gets much smaller than it. For 
these models, the FOTH appears just before the trapping horizons cross, although this is not so clear in the figure. 
For models where the initial BHEH is smaller than the Friedmann cosmological apparent horizon, i.e.. Model E, the 
situation is standard. In terms of r, the FOTH is inside the BHEH, while the POTH is outside it. The BHEH soon 
gets much smaller than the POTH. 

Figure 1^ shows the evolution of the Misner-Sharp mass of the BHEH and trapping horizons. For models where 
the initial BHEH is larger than the Friedmann cosmological apparent horizon, the mass of the BHEH first decreases 
and then increases. As we will see later, the mass of the BHEH decreases if and only if it is in the past trapped 
region. When the black hole gets out of the past trapped region, its mass starts to increase. The mass of the trapping 
horizon with r^y = is only slightly smaller than that of the BHEH, despite the radii being considerably different. 
This is because the density inside the perturbed region is very small. After the BHEH crosses the trapping horizon 
with r^u — 0, which is a POTH at the crossing, the mass of the BHEH soon gets much smaller than the mass of the 
cosmological apparent horizon. 

The change in the BHEH mass may be seen more clearly in Fig.^l which shows the rate of mass increase dmBHEii/dv. 
For models other than E, this is initially negative but it then increases, crosses zero, reaches a positive maximum 
and then decreases. The larger the initial mass ratio of the BHEH to the Friedmann cosmological apparent horizon, 
the more negative the initial mass increase rate is. For Model E, where the initial BHEH is slightly smaller than 
the Friedmann cosmological apparent horizon, the accretion rate starts with a very small positive value, increases to 
a maximum and then decreases. The behaviour of the BHEH mass is clearly explained by Eq. H2.7|l . which can be 
rewritten as a black hole mass equation: 

=const 



dt V dt 



where t = t{u + v) is any time coordinate which depends onu + v. Whether the black hole mass increases or decreases 
depends solely on the sign of r^„, which is the null expansion along the v = const direction. In the usual situation. 
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the BHEH is in a region where r.„ < and its mass monotonically increases. However, when the BHEH is in a 
past trapped region, its mass monotonically decreases because r ^ > 0. Equation (|4.1|) (or H2.7|l ') also explains why 
the mass accretion rate starts very small for Model E. Since the BHEH is inside but very close to the cosmological 
apparent horizon, starts off negative but very close to zero. After some evolution, decreases well below zero 
and the accretion then increases. This suppression of accretion for a PBH as large as a cosmological apparent horizon 
is due purely to general relativistic effects, or relativistic cosmological expansion. Paper I discusses the qualitative 
difference between PBHs whose size is comparable to and much smaller than the cosmological apparent horizon. 

For models with the initial BHEH larger than the Friedmann cosmological apparent horizon, the areas and masses 
of both trapping horizons decrease as v increases before they cross each other. After the crossing, they increase. This 
behaviour is completely consistent with the "second law" for trapping horizons, as formulated by Hayward f^. The 
theorem states, for example, that if the null energy condition holds, then the area and mass of the FOTH (POTH) with 
r^y = do not decrease (increase) along the vector z tangent to the horizon. This vector has the form z = (3dy — adu 
where (3 > 0. 

D. Conversion from a white hole to a black hole 

To understand the spacetime structure for solutions in which the PBH is larger than the cosmological apparent 
horizon, we concentrate on Model F. Figure |H1 gives the detailed numerical results for this model. Figures |HIa) and 
(b) give the evolution of r along the null geodesies with v = const and u — const, respectively, while Figs.|Hl(c) and 
(d) give that of 2m/ r along the null geodesies with v = const and u = const, respectively. It is seen that r continues 
to increase along both the earlier null geodesies with u — const and v = const. This is also consistent with Fig. ^ 
(a). On the initial null ray v — vq — 1, 2m /r decreases monotonically as u increases and crosses unity from above, as 
can be seen in Fig. EJc). Since we have chosen = const for u > Um + Au on the initial null ray v = vq, Eq. H2.6|) 
implies m ~ const here. Since = implies 2m /r = 1, the area radius along the null ray v — vq must monotonically 
increase as u increases. This means that this null ray does not reach r = but another null infinity, different from 
the one which the null rays with u = const reach. 

We can now sec that the standard conformal diagram for PBHs, which is given by Fig. [3 will not apply if the 
BHEH is larger than the cosmological apparent horizon. Rather the conformal diagram is given by Fig.|5| The FOTH 
and POTH intersect at one point in this diagram. This point corresponds to a sphere which is a bifurcating marginal 
surface. It is clear that the spacetime can be interpreted as the conversion of a white hole to a black hole, since a 
PBH larger than the cosmological apparent horizon involves a past trapped region changing to a future trapped one. 
There are two distinct asymptotic regions: the future null infinity for rays with u = const and the one for rays with 
V — const. These two asymptotic regions are associated with two disconnected untrapped regions. It is impossible 
for an observer to go from one untrapped region to the other by traversing the crossing point, the future trapped or 
past trapped region. The spacetime has no regular centre. 

Figure|51is naturally interpreted as an amalgamation of the Kruskal diagram for the extended Schwarzschild solution, 
shown in Fig. 1101 and the standard conformal diagram for a black hole in a flat Friedmann background, shown in 
Fig. 121 It will be recalled that the Kruskal diagram also contains a white hole changing into a black hole and two 
asymptotically flat regions. The analogy of this extension for the conformal diagram shown in Fig. [5] is shown in 
Fig. This contains the big bang singularity, the black hole singularity and two asymptotically flat Friedmann 
regions. There are two BHEHs and two POTHs, these crossing in the black hole region. The regions below the two 
POTHs are past trapped but not precisely equivalent to a white hole since, unlike the situation in Fig. 1101 there are 
no past null infinities and no past event horizons. Since the spacetime has no regular centre, it cannot contain a 
particle horizon as a light cone which emanates from the regular centre just after the big bang. Figure El is obviously 
contained within Fig. llll but the present numerical calculations do not determine what happens in the region v < vq. 

V. DISCUSSION 

If we accept both the conventional conformal diagram for PBHs shown in Fig. |2 and the new one shown in Fig. 
or Fig. 111! we have two distinct PBH causal structures. The conventional diagram has a regular centre r = and 
a single null infinity, while the new one has no regular centre and two null infinities. This raises the issue of what 
determines the causal structure of PBH spacetimes. If we consider standard PBH formation from initial data with a 
regular centre, then the causal structure will be described by the conventional diagram. 

Another issue concerns the transition between the two diagrams. If this transition is governed by parameters which 
can be changed continuously, it would be natural to assume that there is a threshold spacetime which separates them. 
We should note that the BHEH is null, while the cosmological apparent horizon is spacelike in the massless scalar 
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field case (see Fig. O. Hence, if the ubheh is smaller than u = ucph of the cosmological particle horizon (0 in the 
present coordinates), the BHEH can be always outside the cosmological apparent horizon if we take Vq sufficiently 
small. In this case, mBHEH(wo) > "^cah(wo) and hence the causal structure should be given by Fig. or Fig. El On 
the other hand, if ubheh > ^cph, the causal structure is depicted by Fig. |21and the spacetime has a regular centre 
before the BHEH appears. Therefore, the transition spacetime corresponds to a PBH whose event horizon coincides 
with the cosmological particle horizon of the background cosmological solution. The causal structure of this critical 
spacetime is depicted in Fig. 1121 

The discussion of the mass variation of the BHEH can be understood in the more general spherically symmetric 
context from the counterparts of Eqs. (|2.6() and H2.7() 0,0. In our notation, these equations become 

87rr^ 
87rr^ 

The combination of the two terms in parentheses on the right-hand side of Eq. (|5.2|l . related to null expansions, 
determines the time variation of the BHEH mass. In the case of a massless scalar field, the situation is simplified 
because Tuv = and Tw — > 0. We therefore conclude that the black hole accretion vanishes when the event 

horizon coincides with the POTH with r „ = and that it becomes negative when the event horizon is in a past 
trapped region. For general matter fields, it can be proved that the BHEH mass is non-decreasing if r > 0, < 
and the dominant energy condition hold on the event horizon (c.f. Proposition 5 of 7J). However, if these assumptions 
are not satisfied, the mass variation of the BHEH is non-trivial. If the dominant energy condition holds, then the 
conditions Tuv > and T^v > are necessarily satisfied. In the general case with T^v > 0, the BHEH mass still 
increases even when it is on a POTH with r,„ = 0. When the BHEH is in a past trapped region, the first and second 
terms in parentheses on the right-hand side of Eq. H5.2|l make positive and negative contributions, respectively, to the 
black hole mass change. Even when the BHEH is in a past trapped region, its mass can still increase if the first term 
is dominant. This suggests that the BHEH has to be well outside the POTH in order to have suppressed accretion 
or mass decrease for general matter fields satisfying the dominant energy conditions. On the other hand, there could 
still be a lot of accretion for some other kinds of matter field. 

We have shown that the mass of the BHEH can decrease even if all possible energy and causality conditions 
hold. This suggests that the usual concept of event horizons is not appropriate for a proper discussion of black 
hole thermodynamics. This is why Hayward 0,IM0j0| introduced the concept of the FOTH (POTH) as a useful 
generalisation of the idea of a black hole (white hole) event horizon. As we saw in Section IIIBI the idea is that all 
definitions should be quasi-local and this is useful in proving black hole properties. For example, when the black hole 
(white hole) horizon is defined as a FOTH (POTH), one can prove the monotonicity of the area and mass in the 
spherically symmetric situation. From this point of view. Fig. El describes the conversion of a white hole to a black 
hole by definition. 

Although the analysis in this paper has assumed spherical symmetry, it is likely that similar considerations apply 
more generally. If we consider slightly nonspherical PBH spacetimes, it is trivial to show that the strange properties 
discussed in this paper still hold. However, when the system is very far from spherical symmetry, we need to re-analyse 
the Einstein equations. Although the conditions for suppressed accretion or mass reduction may change, we believe 
the appearance of these properties is robust. 

Since a spacetime which contains a PBH larger than the cosmological apparent horizon cannot have a regular 
centre, such PBHs could not form through classical processes from a Friedmann background. This implies that a 
PBH formed through classical processes always increases its mass and soon gets much smaller than the cosmological 
horizon even if it is comparable with the cosmological scale at its formation. On the other hand, if we take quantum 
processes or inflation preceding to the massless scalar field dominated stage into account, the formation of such PBHs 
may be allowed. Moreover, there is a possibility that the universe initially contained a past trapped region which was 
converted into a future trapped region, as indicated in Fig. ^2 This kind of PBH might be described as "ab initio" 
or "eternal" . 



VI. CONCLUSION 

We have investigated PBHs in a flat Friedmann universe with a massless scalar field which are larger than the 
cosmological apparent horizon. Through fully general relativistic numerical calculations, we have found that there 
are two trapping horizons which cross. This corresponds to the conversion of a white hole into a black hole and to an 
initial reduction in the mass of the black hole event horizon. We have seen that these unusual features are peculiar 
to PBHs larger than the cosmological apparent horizon. 
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The present analysis is confined to the spherically symmetric situation with a massless scalar field, so it is clearly 
important to examine how sensitive these peculiar features are to the assumed symmetry and type of matter field. In 
particular, it is interesting to study whether the condition for the black hole mass reduction applies more generally. If 
so, this will be important for black hole physics and black hole thermodynamics. It is also interesting to study under 
what circumstances such a PBH can arise in the early universe, especially in the context of inflationary cosmology, 
and how this relates to the separate universe condition. To make such an extension numerically or analytically, it 
would be very important to adopt the appropriate definition of the cosmological and black hole horizons and the 
trapping horizon framework should be suitable for the case in which the black hole is as large as the cosmological 
horizon. We will study these issues in a future paper. 

The formation of PBHs arising from the evolution of an effectively massless scalar field during inflation and its 
implication for the formation of galaxies and supermassive black holes in galactic nuclei has been recently discussed 
in Rcf. 13] . In this scenario, super-horizon PBHs may naturally arise due to the closed domain walls of size exceeding 
the cosmological horizon. 

Finally, it should be noted that the conversion of a wormholc to a black hole in a first order phase transition was 
studied in Ref. jl^] in the context of bubble nucleation. They showed that if a wormhole is created in the transition 
from a false to a true vacuum, it necessarily collapses to a black hole. Since our system only contains a scalar field 
with no potential and therefore no phase transition, the situation is rather different. However, it is quite interesting 
that similar phenomena are seen in different systems. 
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Initial Singularity 



FIG. 1: The conformal diagram for a flat Priedmann spacetime with a massless scalar fleld. The cosmological apparent horizon, 
which is a past outer trapping horizon (POTH), is spacehke and outside the cosmological particle horizon (CPH). 



BH Singularity 




Initial Singularity 



FIG. 2: The conformal diagram for a PBH smaller than the Friedmann cosmological particle horizon. 
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FIG. 3: Schematic figures for the initial data corresponding to PBHs (a) smaller and (b) larger than the cosmological apparent 
horizon (CAH). 
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FIG. 4: Positions in (m, w) plane of black hole event horizon (BHEH) and trapping horizons with r = and r,u = 0, plotted 
with solid, dashed and dotted-dashed lines, respectively, The future and past outer trapping horizons are labelled as "FOTH" 
and "POTH", respectively. The signs of (r,„) and (r,„) are also shown. A region is future trapped if it has (— , — ), past trapped 
if it has (+, +) and untrapped if it has (+, — ) or (— , +). 




FIG. 5: The area radii of the black hole event horizon (BHEH) and trapping horizons with r,„ = and r^u = 0, plotted with 
solid, dotted and dashed lines, respectively. We can see the area radius of the black hole event horizon always increases. 
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FIG. 6: The Misner-Sharp masses for the black hole event horizon (BHEH) and trapping horizons with r,„ = and r,u = 0, 
plotted with solid, dotted and dashed lines, respectively. The mass of the black hole event horizon decreases when it is in a 
past trapped region and increases when it is in a untrapped region with null expansions (+, — ). The mass of the trapping 
horizon with r,i, = is slightly larger than that of the black hole event horizon for each model, although they are almost 
indistinguishable. 
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FIG. 7: Dependence of mass increase rate dmBHEH/du of the black hole event horizon on v for Models E-J. It is negative when 
the black hole event horizon is in the past trapped region but positive when it gets out of the past trapped region. 
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FIG. 8: The evolution of r along (a) v = const and (b) u = const and of 2m/r along (c) v = const and (d) u = const for 
Model F. In (a) and (c), the arrow at m = —0.5 denotes the matching surface, while the arrow at -u ~ 0.400 denotes the black 
hole event horizon. 
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Initial Singularity 



FIG. 9: Conformal diagram showing the possible causal structure for a PBH larger than the Friedmann cosmological apparent 
horizon. There are two distinct future null infinities and a past trapped (white hole) region is converted into a future trapped 
(black hole) region. The region with w < Wm is the usual flat Friedmann spacetime but the region v < vo is not calculated. See 
text for details. 




FIG. 10: Conformal diagram of the maximally extended Schwarzschild spacetime, i.e., the Kruskal diagram. The black hole 
event horizon (BHEH) coincides with the future outer trapping horizon and the white hole event horizon (WHEH) coincides 
with the past outer trapping horizon. 
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Initial Singularity 



FIG. 11: Conformal diagram of a possible causal structure of the maximally extended spacetime for a PBH larger than the 
Fricdmann cosmological apparent horizon. There are two distinct future null infinities, corresponding to the conversion of a 
past trapped (white hole) region into a future trapped (black hole) region. 



BH Singularity 




Initial Singularity 



FIG. 12: Conformal diagram of the critical PBH spacetime, where the black hole event horizon (BHEH) coincides with the 
cosmological particle horizon (CPH). 



